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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 
RESEARCH MEMORANDUM 

ANALYSIS OF WING LOADS MEASURED ON A FLEXIBLE SWEPT-WING 
JET BOMBER DURING PUSH-PULL MANEUVERS 
By Patrick A. Gainer and Paul W. Harper 

SUMMARY 


Loads on the wing of a Boeing B-47A medium bomber were measured 
by means of strain-gage instrumentation during 23 gradual symmetrical 
push-pull maneuvers performed at altitudes from 15*000 to 30*000 feet* 
Mach numbers from 0.4-7 to 0.8l* and gross weights from 105*500 to 
114*500 pounds . These measurements* together with supplementary measure- 
ments of load factor and other flight parameters* were reduced in order 
to obtain time histories of aerodynamic shear* bending moment* and torque 
at each of four spanwise stations on the left wing and at one station on 
the right wing. These time histories were analyzed in order to deter- 
mine additional and basic load components for the purpose of checking 
currently available methods of predicting loadings on flexible wings. 

From the comparison of experimental and theoretical results* it is 
concluded that presently available methods are capable of predicting 
additional and basic loads accurately within the range of flight condi- 
tions covered by the tests. 

INTRODUCTION 


The quantitative checking* by load measurements in flight* of the 
loads predicted by available calculation procedures continues to be of 
importance* particularly for the case where large interference and flexi- 
bility effects exist. The National Advisory Committee for Aeronautics 
has recently completed a program in which deflections and loads over 
various parts of the relatively flexible Boeing B-47A airplane were meas- 
ured for the specific purpose of providing such data. 

Some phases of this program have been reported in references 1 to 4. 
In the present report the primary consideration is with that phase of the 
program dealing with the measurements of the shear, bending moment, and 
torque at a number of stations across the wing span and the comparison 
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of these results with those predicted by applicable engineering methods. 
The data for this report were obtained during symmetrical push-pull 
maneuvers performed at a slow rate so that quasi— static conditions pre- 
vailed. The maneuvers covered a Mach number range from 0.47 to 0.8l, an 
altitude range from 15,000 to 30,000 feet, and a gross airplane weight 
range of 105 ? 500 to 114,500 pounds. 


SYMB0IS 


b 

Lp 

mR 

M 

n 

q. 

s 

F s 

Mb 

T 

W 

x 


y 


X' 


wing span, in. 

horizontal-tail aerodynamic load, lb 

lift-curve slope of the rigid airplane without the tail, 
per deg 

Mach number 

normal load factor at airplane center of gravity, g units 
dynamic pressure, lb/sq ft 
wing area, sq ft 

aerodynamic shear force, positive for up load, lb 

aerodynamic bending moment, positive for up load, in- lb 

aerodynamic torque, negative for up load rearward of strain- 
gage reference station, in- lb 

airplane gross weight, lb 

coordinate measured from intersection of wing front spar and 
airplane center line, parallel to center line, positive 
forward, in. 

coordinate measured from intersection of wing front spar and 
center line, perpendicular to center line, positive outboard, 
in. 

streamwise distance from a strain-gage reference station to 
the center of pressure of the additional airload carried 
outboard of the station, in. 


NACA RM L57E28 


3 


y* spanwise distance from strain-gage reference station to the 

center of pressure of the additional airload carried out- 
board of the station, in. 

a,Q wing-root angle of attack at zero lift, measured at center 

line, deg 

a r wing-root angle of attack measured at center line, deg 

A sweep angle of wing quarter-chord line, deg 

0 pitching angular velocity, radians/sec 

0 pitching angular acceleration, radians/sec^ 

Subscripts : 

add due to additional load 

av average 

cl referred to the airplane center line 

0 at zero load factor 

sta referred to a strain-gage station 

APPARATUS AND TESTS 


A standard B-V 7 A medium bomber was used to obtain the test data 
analyzed in this report. The only external modifications consisted of 
the addition of a faired nose boom containing the airspeed head and the 
angle -of -attack and sideslip vanes, and the addition of a housing for 
the optigraph (deflection recording) cameras rearward of the cockpit 
canopy. These additions were not considered to cause any appreciable 
load-distribution changes. Figure 1 is a photograph of the test air- 
plane. A three-view drawing showing the modifications is given in fig- 
ure 2. Pertinent physical characteristics of the airplane are given in 
table I. 


Instrumentation 

NACA standard instruments were used to record airspeed, altitude, 
fuselage angle of attack, rotational velocities and accelerations, linear 
accelerations, and control- surface displacements. The distribution of 
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the normal load factor along the wing span was measured hy means of 
remotely recording accelerometers, located as shown in the plan-form 
view of figure 3. Normal and transverse load factors at the nacelles 
were also recorded. 

Electrical wire resistance strain gages were mounted at four span- 
wise stations on the left wing, at one station on the right wing, on 
the engine nacelle struts, and at the horizontal -tail root. The refer- 
ence points for the wing gage stations were located on the front spar 
at 7.8, 36.2, 59.5, and 82.2 percent of the semispan on the left wing 
and at 7.8 percent on the right wing. The gages for each station were 
calibrated by the methods described in reference 5. Equations were 
obtained from the calibrations to give, in terms of bridge outputs, 
the structural shear at each reference station and the structural 
bending moment and torque about the reference axes at each station. 

The bending -moment axes were parallel to the airplane center line, and 
the torque axis for each station was perpendicular to the corresponding 
bending -moment axis. The bridge outputs were recorded in flight on 
multichannel oscillographs. 

Loads applied to the inboard engine nacelle produced some localized 
twisting of the front spar and thereby affected the responses of the 
strain gages in the wing at 36.2 percent semispan. Errors in the wing- 
plus -nacelle shear, bending moment, and torque measured at this station 
were found to be proportional to the nacelle structural pitching moment. 
Provision was, therefore, made for correcting these errors by using the 
measured nacelle pitching moment. 


Tests 

The maneuvers analyzed for this report (runs 5 to 27 of flight 25 
of the general test program) consisted of 23 push-pull or roller-coaster 
maneuvers performed at Mach numbers up to 0.8l and altitudes from 15 , 000 
to 30,000 feet. The gross airplane weight during these maneuvers varied 
from 105,500 to 114,500 pounds. All fuel was carried in the fuselage. 
All pertinent quantities were recorded continuously during each run. 

The maneuvers used for the analysis covered the Mach number range 
in such a way that the Mach numbers were repeated at each of four alti- 
tudes, except that a Mach number of 0.8l could not be attained at the 
lowest altitude (15,000 feet) because of the 425-knot placard speed. 

The gradualness of the maneuvers used was desirable for minimizing 
the effects of pitching acceleration, and because of this gradualness 
some changes in dynamic pressure and Mach number occurred in each run. 
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Table II gives the flight conditions, including the average, the 
maximum, and the minimum dynamic pressure and Mach number for each of 
the 23 runs analyzed. 


METHOD AND RESULTS 


For the maneuvers listed in table II, time histories were recorded 
of structural shear, bending moment, and torque at the four wing sta- 
tions. Time histories were also recorded of all other pertinent loading 
parameters such as load factors, dynamic pressure, Mach number, pitching 
velocity, and acceleration. 

Time histories of aerodynamic shear, bending moment, or torque at 
any wing strain-gage station were obtained by adding the measured inertia 
component to the measured structural loading component. The structural 
component was determined from the strain-gage -bridge outputs by use of 
the calibration equations. The inertia component was obtained from the 
known dead-weight distribution of the wing and the time histories of the 
measured load-factor distribution. 

The measured loadings at the 36. 2 -per cent -semispan station were to 
be corrected for local strain distortions caused by inboard -nacelle air- 
load and inertia load. As stated previously, the intention was that 
the strain gages on the nacelle strut be used for determining these 
corrections at the 36.2-percent-semispan station. The nacelle strain 
gages failed to operate satisfactorily in flight, and it was therefore 
necessary to base the corrections on the nacelle inertia load alone. 
Neglect of the nacelle airload may have contributed a maximum error of 
about 3 percent to the additional aerodynamic loads at the 36.2-percent 
station . 

Failure of the nacelle strain gages also made it impossible to deter- 
mine separate wing and nacelle aerodynamic loadings; therefore, the 
aerodynamic quantities presented in this report are wing -plus -nacelles 
values . 


Analysis of Time Histories 

In the analysis, the aerodynamic -load time histories were considered 
to consist of basic and additional load components as well as possible 
minor load components caused by pitching velocity and pitching accelera- 
tion. The basic load components are those at zero normal load factor 
for the airplane, zero pitching velocity, and zero pitching acceleration. 
The additional component is the increment in loading per unit change in 
airplane normal load factor alone. For this report the primary objective 
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was to determine the basic and additional components of wing shear, 
bending moment, and torque at each station. 


Provided dynamic pressure and Mach number are constant, the fol- 
lowing equations describe the aerodynamic loadings at each station in 
terms of load factor n, pitching velocity 9, and pitching accelera- 
tion 0: 


s, sta ~ 


(•^s, sta) 


^b , sta “ j^b , sta) 


0 


0 



dF 


dF 


n + 


^ sta 0 + - .^ sta 0 (la) 


60 


x •• 

60 


n + q + d M b ; sta £ 


60 


60 


(lb) 


T 


sta 


( T sta) 


6T, 


sta 


0 


6n 


n + 


6T, 


sta 


60 


e + 


6T, 


sta Q 


60 


(lc) 


The quantities enclosed by brackets in equations (l) are, in the 
order of appearance, the basic and additional components of aerodynamic 
load. Equations (la) to (lc) were # solved by least-squares procedures by 
using the measured values of n, 0, 0, and so forth, to obtain the 

bracketed quantities for each strain -gage ## station and each run. The 
components of loading caused by 0 and 0 were found to be negligible 
for the runs analyzed, and these results are not presented in this report. 


(When the dynamic pressure and Mach number are variable, the coef- 
ficients of n, 0, and 0 and the intercepts ^F g s ^_ a j , ^M-^ , 

and ^ s -^ a ) are ncy k constants. Since each of the runs analyzed con- 

tained some changes in q and M, a trial analysis of one run was made 
by the method described in appendix A to determine whether changes in 
q and M during a run should be taken into account. Equations (la), 
(lb), and (lc) were found to be sufficiently accurate to determine values 
of the coefficients and intercepts corresponding to the average q and 
M of each run . ) 


In addition to the shears obtained from the wing strain-gage meas- 
urements, the total lift of the wing-fuselage -nacelles combination was 
determined from the known airplane weight W, the airplane normal load 
factor n, and the measured tail load Lp. The total load for the wing- 
fuselage -nacelles combination is simply nW - Lip. The aerodynamic shear 
at the airplane center line is one -half the total normal load. The time 
history of the center -line shear for the wing-fuselage-nacelles combina- 
tion was separated into basic and additional components by the method 
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used for shears at other stations , expressed hy the following equation, 
which is similar to equation (l): 


i(nW - L T ) = ifnW - L T ) + 1 d ( nW ~ n + i ^ ' ■ i + 1 ^ ■ ' 

2' 2 \ T /o 2 dn 2 d e 2 de 


or 


+ 


s, cl 


2 x 0 dn 2 de 


( 2 ) 


1 dFq p 7 -| dfnW - Lm ) • -| d^nW - Lm) •• . . 

i Lm + -- S J C > - n + i -A — _ — IZ 0 + A _A _JEZ. 0 (2a) 

o q Aq 2 


The term 


dF 


s ,cl 


dn 


represents the additional component of aerodynamic 


shear at the airplane center line due to the lift of the wing-fuselage- 
nacelles combination. The additional aerodynamic shears, bending moments, 

q C* 7 

and torques at each of the strain-gage stations, when divided by — - 


produce the quantities 


dn 


Shear ratio = 


s , st a ' 


s ’ cl ) add 


ratio 


, sta\ 

Fs > cl /add 


Torque ratio = 



The quantities are the ordinates for the plots of figure Parts (a), 
(b), and (c) of figure 4 give, respectively, the shear, bending moment, 
and torque at each strain-gage station per unit of additional aerodynamic 
shear at the center line of the wing-fuselage -nacelles combination. The 
curves in these plots represent the predictions of these quantities made 
by the method outlined in appendix B. 

Because the theoretical calculations were made in terms of the 
parameter qm^, it was considered appropriate to plot both predicted 
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and experimental results against qm R for purposes of comparison and 

evaluation. Therefore, the abscissas for the experimental ordinates in 
figures A and 5 are the values of the product of the average dynamic 
pressure q and the lift-curve slope m R for the rigid airplane with- 
out the tail corresponding to the average Mach number of each run. The 
values of m R were obtained from the curve of m R against M derived 

from wind-tunnel tests and given in reference 6. 


Spanwise and streamwise distances from the strain-gage reference 
stations to the centers of pressure of the additional airload outboard 
of each station were obtained from the data by dividing additional- 
airload bending-moment and torque ratios by the corresponding shear 
ratios. Thus, 


The coordinates y^ ta and x^ ta are plotted against qm R in fig- 
ure 5. The predicted center -of -pressure coordinates obtained as described 
in appendix B are represented by the curves in figure 5. 

Figure 6 shows the locations of the centers of pressure of the addi- 
tional airload outboard of the 7. 8-percent -semispan station, as they 
appear on a plan -form view of the wing. Further use was made of the 
additional aerodynamic shears obtained by the use of equations (la) and 
(2a) in order to calculate wing deflections due to additional airload. 
These calculations were made for only one run for the purpose of com- 
parison with some of the results of deflection measurements reported in 
reference 2. Run 16 (table II) was found to be nearly equivalent to a 
maneuver (flight A, run 20) of reference 2 except for a slight differ- 
ence in gross airplane weight. The additional shears at each station 
obtained from the analysis of run 16 were adjusted for this weight 
difference . 


Additional -Load Centers of Pressure and Deflections 
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Deflections at each of eight optigraph target stations were calcu- 
lated from these adjusted shears by the use of an influence -coefficient 
matrix derived from the experimental data reported in reference 7 • This 
matrix gives the optigraph target deflections due to concentrated loads 
applied at each of nine stations along the wing quarter-chord locus. 

The wing was divided into nine segments, each with one of the loading 
stations at its midpoint. The change in shear over any one segment was 
considered to he a load concentrated at the midpoint of the segment. 

These changes in shear were estimated by plotting the adjusted shears 
from run l 6 against span station, fairing a curve through these points, 
and reading the shears at the segment boundaries . The average results 
of right and left wings were used at the 7 . 8 -percent-semispan station. 

The deflections thus calculated are those due to additional airload 
only, while those measured in flight include the deflections due to iner- 
tia load. The deflections due to inertia load are given in reference 2. 
These were subtracted from the deflections given for ruin 20 in order to 
obtain deflections due to additional airload. 

The deflections calculated from the measured loads are compared 
with the results of deflection measurements in figure 7 > where both sets 
of deflections are plotted against span station. The missing points in 
the measured deflections are the result of failure of two of the opti- 
graph lights. 


Basic Loadings 

The intercepts of the measured aerodynamic loadings at zero load 
factor, the first term in equations (l), do not necessarily represent 
the true basic shears, bending moments, and torques. This fact is illus- 
trated in figure 8 , which shows the left-wing root -shear intercepts at 
zero load factor plotted against qm R for the different altitudes of 

the test maneuvers. A curve faired through the points at any single 
altitude would be expected to pass through the origin of the plot because 
there can be no aerodynamic load if qm R is zero. Obviously, this con- 
dition is not met in figure 8 ; rather, the points at each altitude appear 
to lie on a nearly straight line, the line for one altitude being approx- 
imately parallel to the lines for the other altitudes. Similar plots of 
the other basic load quantities showed the same kind of zero shifts. 

An equation representing all the basic shears at one station by a 
function of qm R , while allowing a different intercept for each altitude, 

is as follows: 

( F s,sta ) 0 = A ( qm R) + B ( qin R ) 2 + ^ 5,30 + ^3,25 + ^ 3,20 + ^8,15 


(3) 
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The AF S terms in equation ( 3 ) are considered to be the zero 

shifts in the measured load with changes in altitude , and the numerical 
subscript indicates the altitude (in thousands of feet) to which each 
pertains . 

Values of AF S obtained by least -squares procedures were applied 

as corrections to the basic shears. Similar corrections were made to the 
basic bending moments and torques. 

The resulting adjusted basic loads are presented in figure 9- Parts 
(a), (b), and (c) of figure 9 show, respectively, the basic shears, 
bending moments, and torques from equations (l) adjusted by the method 
described by equation 3 for each strain-gage station and each run. The 
abscissas are again the average values of qm^ for each run. The basic 
shears, bending moments, and torques predicted by the method outlined in 
appendix B are given by the curves in figure 9 . 


DISCUSSION 


Additional-Load Quantities 

The curves of figure 4 show that the predicted additional aerody- 
namic loading quantities agree well with those determined from flight- 
test measurements. In each case, the experimental results form a curve 
which is approximately parallel to the predicted curve, though in some 
cases there are noticeable differences in level. These differences in 
level may be due to slight errors in either the flight measurements or 
the wind-tunnel measurements on which the predictions were based. The 
possibility of small measuring errors is indicated by the consistent 
differences between the measurements at 7.8 percent semispan on the right 
and left wings, especially in figure 4(a); however, because the measured 
right and left shear ratios are both below the predicted curve, it is 
possible that the wing outboard of the fuselage carried somewhat less 
than the predicted fraction of the total load. From the wind-tunnel 
pressure-distribution tests, the exact fraction of load carried by the 
fuselage was indeterminate. The load-distribution curve from that source 
was, therefore, faired smoothly into the center line to obtain the pre- 
dicted results at zero qmp. 

Small discrepancies in level at zero qm^ apparently have little 

effect on the predicted changes in the additional -load shear, bending 
moment, and torque ratios, because the predicted results approximately 
parallel the experimental results, even at the highest values of qm-p. 
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The small scatter in the flight results is apparently random. No 
correlation could be found between changes in airplane weight and the 
displacement of the points from an average curve. This result was to 
be expected because of the relatively small changes in airplane weight. 

The quantities given in figure 4 are considered to be primary quan- 
tities. From these quantities the center -of -pres sure coordinates in 
figure 5 were obtained. The discrepancies in figure 5 between predicted 
and experimental results are associated with discrepancies previously 
noted in connection with figure 4. The greatest of these discrepancies 
are at 7.8 percent semispan, where there are large differences between 
experimental results for right and left wings and at the 36 - 2 -percent 
station. All major differences between predicted and experimental center - 
of -pressure coordinates are seen to be the result of relatively smaller 
differences in the quantities of figure 4. 

The center-of -pres sure results for the wing-root station may be 
taken as representative of the changes in additional-load center of 
pressure of the entire wing -body -nacelle configuration caused by flexi- 
bility. Figure 6 shows the locations of the centers of pressure of the 
additional airload outboard of the 7 . 8 -percent-semispan stations, as 
they appear on a plan-form view of the wing. Only that section of the 
wing is shown which contains the centers of pressure and M.A.C.^, the 

mean aerodynamic chord of the area outboard of 7-8 percent of b/2. The 
quarter-chord locus is given also. Experimental results for both right 
and left wing have been plotted together for easy comparison, and the 
predicted results for the flight range of qm^ are given by the curve. 

The direction of movement is inboard and forward with increasing qm.p, 

as designated by the arrow. Because of the effect of the inboard engine 
nacelles, both experimental and theoretical centers of pressure lie on 
loci somewhat forward of the quarter-chord locus. Had the effects of 
the nacelles been neglected in predicting the loadings, the predicted 
centers of pressure would have fallen on the quarter -chord locus. 

Although the initial location of the center of pressure (at the 
lowest value of qmp) is in doubt because of the discrepancies between 

experiment and prediction and also between the two sets of experimental 
results in figure 6 , it is fairly well established that the changes in 
additional -load center of pressure were accurately predicted and that 
the effects of nacelle loads are not negligible for this airplane. The 
fact that the nacelle effects were strong enough to give a center-of- 
pressure location noticeably forward of that which would be predicted 
for the wing alone indicates that the nacelles should also have an 
appreciable effect on the changes in center of pressure due to flexi- 
bility. For the change in center -of -pressure location with change in 
qm R , theoretical calculations show that, over the qm R range from 0 to 
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50 , the movement of the center of pressure for the wing alone would he 
about 18 percent greater than that calculated for the wing plus nacelles. 
(A line representing the predicted center-of -pressure locus for the wing 
alone in figure 6 would fall on the wing quarter -chord locus and would 
be approximately 18 percent longer than the line shown.) The combined 
effects of the inboard nacelle are, therefore, to move the center-of - 
pressure locus forward of the wing quarter-chord line and to decrease 
the amount of center-of -pressure movement over a given qm^ range. 


Correlation With Deflection Measurements 

The comparison of wing deflections obtained directly from flight 
measurements with those calculated from measured loads shows a high 
degree of correlation between the two different types of flight-test 
measurements (fig. 7)- Although it is not shown in figure 7 , the sub- 
stitution of predicted shears for measured shears in the deflection 
calculations yielded almost identical results. This agreement of results 
is an indication that the differences between experimental and predicted 
shear ratios in figure 4(a) are not significant where wing-deflection 
calculations are concerned. 


Basic-Load Quantities 

The agreement which was obtained between predicted and measured 
basic shears, bending moments, and torques in figure 9 is fairly good. 

It should be stated that strain gages are not entirely suited to the 
measurement of basic loads because of their response to strains produced 
by changes in temperature as well as to those produced by actual loads. 
The strain gages near the nacelles may be subject to the effects of 
changes in engine temperature, and these changes are nearly proportional 
to qm R . These effects cannot be detected by the method of zero-shift 

correction described in the section entitled "Method and Results"; fur- 
thermore, the theoretical calculation of basic loads is dependent upon 
a number of wind-tunnel measurements of relatively small pitching -moment 
coefficients. Errors in the wind-tunnel measurements of wing and nacelle 
pitching moments would cause relatively larger errors at the outboard 
stations than at the root station. The root -station basic loads are 
mainly affected by the zero-lift pitching -moment coefficient for the 
entire airplane without the tail. 

The maximum basic shears occurring at the root station are about 
7 percent of the maximum structural shear for which that station was 
designed. While the basic loads are not negligible from the structual 
design standpoint, an error as great as 20 percent in their determina- 
tion could probably be tolerated. The actual discrepancies between 
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predicted and experimental root -station basic loadings are much less 
than this . 

The greatest differences, percentagewise, occur at the two stations 
near the engine nacelles . This might be expected from the previous 
statements about strain-gage response to temperature gradients and the 
dependence of the calculated loads at the outboard stations on the meas- 
urement of relatively small quantities by wind-tunnel testing. 


CONCLUSIONS 


Aerodynamic shears, bending moments, and torques measured on the 
wing of a B-47A airplane during quasi -static symmetrical maneuvers have 
been analyzed in order to obtain basic and additional components. The 
values of these components were then compared with values of the same 
quantities predicted by the theoretical extension of low-speed wind- 
tunnel measurements. As a result of the analysis and comparison, the 
following conclusions are drawn: 

1. Additional-load quantities, including centers of pressure, can 
be adequately predicted. Changes in these quantities with increasing 
values of the parameter qm^ can be accurately determined by the theo- 
retical methods available. 

2. Basic shears, bending moments, and torques can be adequately 
predicted, especially near the root. 

3- A comparison of deflections calculated from measured loads with 
deflections measured in flight shows good correlation between the two 
different types of measurements. Both sets of experimental deflections 
agree well with theoretical deflections. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., May 9, 1957. 
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APPENDIX A 

EFFECTS OF CHANGING DYNAMIC PRESSURE AND 
MA.CH NUMBER ON AMLYSIS 


In order to determine the effect on the coefficients of changes in 
dynamic pressure and Mach number during a run, a single maneuver was 
analyzed by two methods. One of the methods was based on the assumption 
that the additional- and basic-load quantities for the average dynamic 
pressure and Mach number of a run would be obtained if changes in q 
and M within the run were neglected. The other method assumed linear 
variations of basic and additional quantities with changes in the param- 
eter qm R . The deviations from the average value of qm R were used to 

establish the average basic and additional quantities for the run. 

Run 26 was used for this purpose. (See table II for variations in q 
and M . ) 

The bending moment at the root station was considered to be most 
likely to show the effects of q and M changes . For the first method 
used, the equation for bending moment in^terms of load factor n, pitching 
velocity 0, and pitching acceleration 0 was written as 


Mb. 


sta 


(^b, sta) 


^ av L 


^b , sta 


dn 


n + 


d % sta e + ^sta q (A1) 


av 


<30 


de 


The terms in brackets are of primary interest. ^ M b,sta^Q 

the basic bending moment when all other quantities are zero, and 
d^b,sta 


dn 


is the additional aerodynamic bending moment. 


av 


is 


av 


Because both bracketed quantities in equation (Al) would be expected 
to vary with q and M as the product of q and the lift -curve slope 
m R of the rigid airplane without the tail, the corresponding equation 

for the second method was written as 


(Mb, sta) 0 



+ K 0 A ( cim R) 


d^b , sta 
dn 


fd M 


b, sta ] 


dn 


av 


K n A(qm R ) 


(A2) 


(A3) 
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where 


A(qm R ) = qm R - (qm R ) 


av 


The quantities Kq and K n are constants if a linear variation of 
dM b -ta 

q and — ^ ■ - - with qm R is assumed over the range of qmp for 
a single run. Substituting equations (A2) and (A3) into (Al) gives 

1 / dM b sta\ 

M b,sta. * («b,sta) 0 av + *0 A (*”e) * ) ” + “n A h*)“ * 

— — l cx v \ / a V 


^*b , sta • , sta 

60 + 60 


(A4) 


For the purposes of the analysis represented by equation (A4), the 
assumption that changes in m R are given by the Glauert factor for 

swept wings is sufficiently accurate. If the lift-curve slope at zero 
Mach number is denoted by ( m Rj M lift -curve slope at any Mach 


number is 


m R = 


Nm=< 


o 


\[T- 


P P 
M cos^A 


and 


qi^ = 


Mm 




M^cos^A 


It is not necessary that m R ^ be known since it can be included in the 
unknown coefficient of A^qm^j in equation (A4) . The quantity 


( mR ^M=0 


\fl - M 2 


p 

cos A 


is substituted for A^qm R j and gives 
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/m \ 

+ (n 1 \ A ^ + 

d ^b , sta 

[v d , sta j qJ 

„ orWjbo ./ 5 ' 2 

av V 1 - M COS^A 

dn 


n + 


MM 


M=0 


V 1- M 2 cos 2 A 


n + ^ . sta 0 + dM t) , sta £ 

de de 


(A5) 


Equations (Al) and (A5) were each solved by least -squares methods by 
using 25 points from the time histories of bending moment, load factor, 
q_, M, 0, and 0 for rim 2 6. The basic and additional components of 
bending moment from each equation and the probable errors of fit are 
given in the following table for comparison: 



^b,sta 


( M b,sta) 0 _ 

) 

Probable error, 

Equation 

dn 

in-lb /g 


in-lb 

av 

in-lb 

(Al) 

12.16 x 10 6 


0.13 x 10 d 

O.O65 x 10 6 

(A5) 

12.21 x 10 6 


.12 X 10 6 

.065 X 10 6 


Because the results from equation (Al) are nearly the seme as those 
from equation (A5) , the basic and additional-load components corresponding 
to the average qm R of any run of these tests could be determined accu- 
rately, and with a considerable saving in time, by neglecting changes in 
q and M during a run. 

The coefficients for 0 and 0 were indeterminate because of the 
very small changes in these quantities. 
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APPENDIX B 

THEORETICAL ADDITIONAL AND BASIC LOADS 


The additional aerodynamic wing-load distribution was calculated 
by considering separately the effects of angle -of -attack and load -fact or 
changes. The matrix method given in reference 8 was used in combination 
with the superposition method of solution given in reference 9- This 
method of solution was modified to incorporate least -squares methods of 
curve fitting in order to define the elastic twist distributions in 
terms of span station. 

The required aerodynamic and structural influence -coefficient 
matrices were calculated for nine stations along the semispan. The 
spanwise distribution of section-lift -curve slopes for the rigid wing 
was determined from low-speed wind-tunnel pres sure -distribution tests 
(ref. 10) by the method of reference 8. 

The equilibrium load distributions resulting from a unit symmetrical 
change in wing root angle of attack were computed for a number of values 
of the product of the dynamic pressure q and the lift -curve slope m-p 

of the rigid airplane without the tail while the load factor was held 
constant. Since both fuselage overvelocity effects and nacelle lift and 
pitching -moment effects are proportional to angle of attack (see ref. 8), 
these effects were superposed on the uniform unit angle of attack. Over- 
velocity effects were estimated from the equations in reference 8, and 
the nacelle lift and pitching moment were estimated from wind-tunnel 
tests. Values of qm R were chosen which would approximately cover the 

range of the flight tests and which would permit defining a smooth curve. 

The equilibrium airload distributions due to wing inertia twist per 
unit airplane load factor were similarly calculated while the root angle 
of attack was held constant. 

The results of the two preceding steps were, for each value of qm R , 

the partial derivative of aerodynamic load with respect to root angle of 
attack and the partial derivative of aerodynamic load with respect to 
load factor, for each of the nine span stations. Successive integrations 
of these load distributions gave the partial derivatives of shear and 
bending moment with respect to angle of attack and load factor. Torques 
were obtained by integration along the streamwise projection of the wing 
quarter-chord locus . Concentrated nacelle loads were added at the proper 
stations in each case. 
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The derivative of the aerodynamic shear, bending moment, or torque 
at a particular wing station with respect to airplane normal -load factor 
is related to the partial derivatives of the particular loading with 
respect to angle of attack and load factor and to the total derivative 
of root angle of attack with respect to load factor. The following 
equation, which applies to shear, is an example of this relationship. 

< ^ F ’s,sta _ dF s,sta ^r + / dF s,sta 
dn da^. dn y dn 


(Bl) 


inertia 


The term 


dF 


s^sta I in equation (Bl) is the airload caused 

inertia 


dn 


by wing twisting due to inertia load and must not be confused with the 
inertia load itself. 

At the airplane center line, the derivative of shear with respect 
to load factor is known to be i/w - ^T^ . At the center line, therefore, 


dn y 

1 L _ ^t\ dF s,cl ^r ( dF s,c l \ 

2\ dn J dc^ dn ( dn jj 


(B2) 


inertia 


rinr. 

Equation (B2) was used to determine for each value of qm R 

for which values of the partial derivatives were calculated and for sev- 

dL„ 

eral values of W 1 which covered the range of the flight-test val- 

dn 

values of thus calculated were inserted in equation (Bl) 


ues . The 


dn 


in order to obtain the additional shear, bending moment, and torque at 

dLn, 

each station for each value of qm R and for each value of W - ^ . 

These shears, bending moments, and torques were nondimensionalized with 

if dL T \ 

respect to the additional center -line shear -IW - j, just as were 

the experimental results . The changes in these loading ratios caused 

dL™ 

by changes in W were found to be negligible, even for a 20- 

dn 

dLrp 

percent change in W The theoretical loading ratios corresponding 

dn 
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to W ^ = 110 , 000 , approximately the average flight-test value, were 

dn 

used for comparison with the flight-test results . 

The predicted additional -load shear, bending -moment, and torque 
ratios were plotted against qmp and appear in figure 4 in comparison 

with the flight-test results. 


At zero load factor for the airplane, the sum of the normal air- 
loads on the airplane must he zero. Simultaneously, the sum of the 
pitching moments about any pitch axis must also be zero. 

At the airplane center line, the shears due to the aerodynamic 
forces on one-half of the airplane at zero load factor are: 


1. Wing shear due to root angle of attack, 


dF 


s , cl 




2. Wing shear due to twisting of the wing by wing and nacelle aero- 
dynamic pitching moments, ^F g ^ 


'cl 


5 . Wing shear due to fuselage overvelocity, described in reference 8 
( F s,over ) d 


4. Horizontal -tail aerodynamic load, i Lp q 

5. Nacelle airload, (F s ,na C ) ( 


'cl 


Equating the sum of the vertical forces to zero gives the following 
equation: 


" a °) + ( F s, 0 cZ + ( F s,over ) c z + ( F s,nac) cl + 


0 

(B3) 


Each of the forces in equation (B3) may exert a pitching moment 
about an axis through the quarter-chord point of the mean aerodynamic 
chord. In addition, there will be a pitching moment about this axis due 
to the aerodynamic pitching -moment coefficient q of the airplane 

without the tail. Equating the sum of the pitching moments to zero gives 
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dF 


s, cl 




(“r - °o) + ( F =,t) 


c! + 2 ^ + C M,01S f = 0 (*) 


where Ax' is the distance measured from 0.25 mean aerodynamic chord 
streamwise to the center of pressure of the load indicated in the 
subscript. 

In order to compute the components of equations (B3) and (B4), the 
pitching -moment coefficients of the wing, nacelles, and airplane with- 
out the tail were estimated from the wind-tunnel -test results given in 
reference 6. Each of these coefficients was assumed to vary with Mach 
number as mo does . 


Using the aerodynamic and structural influence-coefficient matrices 
and methods of solution which were used in the additional-load calcula- 

dF, 


tions, distributions of F s ^, 


— , and F 


s , over 


were calculated for 


each of a number of values of qm^. Bending -moment and torque distri- 
butions and the necessary Ax' distances were also calculated. The 
pitching moments due to fuselage overvelocity and engine nacelles do 
not appear in equation (B4) . These quantities are included in the Cjj q 

term as measured in the wind-tunnel tests. 


In both equations (B3) and (B4), only Lip ^ q and a r - cx,q are 

unknown. These quantities were obtained by simultaneous solution of 
equations (B3) and (B4) for each value of qmp. 

At any strain -gage station, the basic shear is 



^s,sta 

doj. 


(Or - Oq) + ( F s,over) sta + 


( F s,nac) 


sta 

(B5) 


All quantities on the right-hand side of equation (B5) are known as 
a result of the preceding computations . 


The bending moment and torque due to load were computed from their 
component parts by equations similar to equation (B5) . The torque due 
to load was added to the torque due to wing and nacelle aerodynamic 
pitching moments in order to obtain the basic torque for each station 
and each value of qmp. 
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TABLE I 

WING DIMENSIONS AND CHARACTERISTICS OF THE B-47A AIRPLANE 


Span, ft Il6.0 

Area, sq ft 1,428.0 

Aspect ratio 9*43 

Taper ratio 0.42 

Thickness ratio (constant) 0.12 

Mean aerodynamic chord, in 155*9 

Sweep of 0.25 -per cent -chord line, deg 35*0 

Root chord, in 208.0 

Tip chord, in 87 . 0 

Airfoil section BAC 145 

Incidence (root and tip), deg 2.75 

Dihedral, deg 0 


TABLE II 


FLIGHT CONDITIONS 


[Center of gravity at 22.5 ± 1 percent mean aerodynamic chord for all runs 


Flight 25 , 
run 

Pressure 

altitude, 

ft 

M 

av 

^max 

“min 

^av' 

lb/sq ft 

^max’ 
lb/sq ft 

‘Wn' 

lb/sq ft 

w, 

lb 

Wav 

5 

50,000 

0.81 

0.823 

0.807 

293 

302 

285 

114,500 

32.8 

6 

30,000 

.76 

.76k 

.760 

258 

259 

255 

114,200 

27-9 

7 

30,000 

.73 

.745 

• 715 

230 

238 

216 

113,800 

24.4 

8 

30,000 

.67 

.683 

.654 

197 

206 

187 

113,400 

20.1 

9 

30,000 

.61 

.626 

.602 

162 

172 

155 

113,100 

16.0 

10 

30,000 

.56 

.581 

• 545 

141 

151 

131 

112,700 

13.5 

11 

25,000 

.81 

.817 

.807 

365 

372 

357 

112,500 

40.9 

12 

25,000 

.77 

.780 

• 755 

335 

349 

318 

112,200 

36.4 

13 

25,000 

.71 

.715 

.699 

278 

289 

268 

111,800 

29.1 

14 

25,000 

.66 

.675 

.647 

244 

259 

230 

111,300 

24.8 

15 

25,000 

.62 

.645 

.596 

215 

237 

193 

110,900 

21.4 

16 

25,000 

• 57 

.587 

.544 

178 

196 

161 

110,600 

17.2 

17 

20,000 

.81 

.818 

.805 

445 

454 

434 

110,100 

49.8 

18 

20,000 

.76 

.760 

• 750 

387 

398 

379 

109,800 

41.8 

19 

20,000 

• 72 

.780 

• 700 

360 

383 

336 

109,400 

37-9 

20 

20,000 

.66 

.670 

.655 

294 

311 

278 

109,100 

29-9 

21 

20,000 

• 57 

• 592 

.540 

225 

251 

200 

108,600 

21.8 

22 

20,000 

.47 

.490 

.428 

145 

168 

123 

107,700 

13.4 

23 

15,000 

• 77 

• 773 

.752 

490 

507 

470 

107,300 

52.9 

2k 

15,000 

• 71 

.716 

.692 

425 

445 

408 

107,000 

44.4 

25 

15,000 

.66 

.669 

.644 

374 

394 

35 ^ 

106,700 

37-8 

26 

15,000 

.56 

.582 

.542 

270 

294 

247 

105,900 

26.1 

27 

15,000 

.1+7 

.494 

.425 

180 

210 

148 

105,500 

16.6 


ro 
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Figure 1.- Test airplane. L-86692 
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Figure 2.- Three vievs of test airplane. 
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Figure 3.- Strain-gage and accelerometer locations. 
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(a) Shear ratios. 

Figure 4.- Comparison of experimental and predicted load ratios. 


Bending moment at each strain-gage station per unit additional aerodynamic 
center-line shear, f ^ ' s ‘ ia ') > (in -lb)/lb 
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(b) Bending-moment ratios. 
Figure 4.- Continued. 
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(c) Torque ratios. 
Figure 4.- Concluded. 
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Figure 5.- Coordinates of centers of pressure of additional aerodynamic 
load against qmp in comparison with predicted coordinates . 
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(b) Strearawise coordinates. 


Figure 5*- Concluded. 
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Figure 6.- Experimental and theoretical locations of additional-airload 
center of pressure with respect to quarter-chord locus and mean aero- 
dynamic chord for area outboard of 7.8 percent semispan. 
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Figure 7 .- Wing deflections calculated from measured additional airload shears compared with 
measured deflections from reference 2. Measured rear-spar deflections not available at 
y = 212 and 68 l. 
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Figure 8.- Experimental basic shears from left-wing root measurements 

illustrating zero shifts. 
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(a) Basic shears. 

Figure 9-- Adjusted basic loads at each strain-gage station in compari- 
son with predicted basic loads. 
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(b) Basic bending moments. 


Figure 9.- Continued. 
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(c) Basic torques. 


Figure 9* _ Concluded. 
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